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Abstract 

We consider the quantum effective action of Dirac fermions on four dimen- 
sional flat Euclidean space coupled to external vector- and axial Yang-Mills 
fields, i.e., the logarithm of the (regularized) determinant of a Dirac operator 
on flat R 4 twisted by generalized Yang-Mills fields. According to physics folk- 
lore, the logarithmic divergent part of this effective action in the pure vector 
case is proportional to the Yang-Mills action. We present an explicit compu- 
tation proving this fact, generalized to the chiral case. We use an efficient 
computation method for quantum effective actions which is based on calcula- 
tion rules for pseudo-differential operators and which yields an expansion of the 
logarithm of Dirac operators in local and quasi-gauge invariant polynomials of 
decreasing scaling dimension. 

MSC-class: 81T13; 58J42; 35S99 

1 Introduction 



Determinants of differential operators arise as (exponentials of) effective actions in 
quantum field theory. The precise definition and investigation of such objects is an 
interesting and challenging mathematical problem which has lead to an active and 
fruitful interplay between mathematics and physics. 

In this paper we compute the logarithmic divergent part, S\ og (A), of the logarithm 
of the regularized determinant for Dirac operators describing Dirac fermions cou- 
pled to a generalized Yang-Mills field A on four dimensional spacetime. For sim- 
plicity we assume spacetime to be flat M 4 with Euclidean signature and the natural 
spin structure. The Yang-Mills fields we consider contain, besides the vector part V, 
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also a chiral (axial) part C (for precise definitions see Eq. (0) ff. below); we write 
A = (V, C). Our definition of S\ og is motivated by physical considerations and will be 
explained further below. To indicate the mathematical significance of our calculation, 
we note that S\ og (A) is (essentially) the noncommutative residue |Wo|| of the loga- 
rithm of Da (see Eq. (|]) for the precise statement). A main motivation for this work 
is to present a computation method for effective fermion actions which at the same 
time is mathematically rigorous, close to standard Feynman diagram computations 
in quantum field theory (see, e.g., |Wej| ), and simple to use. We believe that this 



method is a useful alternative to other methods like the (^-function regularizations or 
the heat kernel expansions (see, e.g., |BGV|| ). We therefore made some effort to 



present this method in a self-contained way, in the hope that this is useful also for 
readers who are mainly interested in learning how to compute effective actions. 

We now discuss our computation method (parts of this method were used pre- 
viously by us in JC], |L1V1| |). We regard the Dirac operator as a PSDO (pseudo- 
differential operator) on a Hilbert space of square-integrable functions on M 4 . Our 
starting point is the following definition for the regularized effective fermion action, 

S A (A) :=Tr A flog (°±p!l) - log ^ + 



Ao / V. Ao 

where m is a real parameter which has the physical interpretation of a fermion mass, 
and A is a positive regularization parameter which we call UV (ultra-violet) cutoff'. 
The role of the non-zero and complex parameter A is two-fold. Firstly, it makes the 
argument of the logarithm dimensionless, and secondly, setting A = |A |/(1 + i0 + ) 
avoids possible ambiguities due to the branch cuts of the logarithm which otherwise 
can arise.[] This definition above has three ingredients. Firstly, a definition of the log 
of an operator a as an integral of the resolvent of a. Secondly, some basic facts about 
PSDO which imply a simple and powerful formula for the symbol of the resolvent of 
the Dirac operator D^. And thirdly, a definition of a regularized Hilbert space trace 
TrA (where removing the regularization corresponds to the limit A — > oo). Combining 
these ingredients we obtain an expansion of Sa(A) in local and quasi-gauge invariant 
polynomials of decreasing scaling dimension. We find 

S A (A) = A 2 S^(A) + log (-^j S log (A) + S<®(A) + C^A" 1 ), (2) 

and this provides our definition of S\ og (A). Our results for S\ og (A) and S^(A) will 
be presented in the next Section. We shall also demonstrate on our way that S\ og (A) 

1 Of course, all results must be independent of |Aq|, and this is a useful check. 
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is proportional to the noncommutative residue [|Wo|| of the logarithm of the Dirac 
operator D^, 

v( , H4Res ( log (^)_ log (2^)). (3) 

The logarithm of the regularized trace of the determinant of the Dirac operator can 
then be defined as 

where TR is the renormalized trace which we will define, and we will provide all 
mathematical tools necessary for computing S^°'(A) explicitly. 

We note that our computation method is closely related to methods which have 
been used in the physics literature for a long time (see, e.g., ||dW| , |Z| , |Wej| ). The 



regularization we use is simple and close to how regularizations are often done in 
Feynman diagram computations, i.e., by introducing a sharp UV cutoff (see Eq. 
(|32l)). We believe, however, that we can offer some improvements in detail which 
make computations easier, more transparent in structure, but nevertheless such that 
each step can be easily justified with mathematical rigor. 

We now discuss some motivation for our computation from a quantum field theory 
point of view. As was known already to Schwinger for the Abelian case, the effective 
action of fermions coupled to a Yang-Mills field A = V (i.e., C = 0) contains a 
logarithmic divergence, log(A/m) S\ og (A), and Si og (.<4) (for C = 0) is proportional to 
the usual Yang-Mills action 



Ym{A) = —j / (i -r u: J-,,,, 

(see, e.g., |Q, Eq. (12.123) where 1/e corresponds to log(A/m)). This is important 
since it implies that a change in the cutoff in the gauge theory, A — > A', leads to a 
finite change of the effective fermion action which can be absorbed by changing the 
Yang-Mills coupled constant, g~ 2 — > (g')~ 2 = g~ 2 + const. log(A'/A). The logarithmic 
dependence of the Yang-Mills coupling constant on the UV cutoff is remarkable and 
distinguishes four spacetime dimensions from all others. 

Our computation is closely related to more recent ideas which have lead to a 
deeper geometric understanding of the standard model of elementary particle physics 
(including Higgs sector). This approach is based on Connes' NCG (noncommutative 



geometry; textbooks on this subject are, e.g., |GVF|| ). One important ingredient 
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of this approach is to define a generalized Dirac operator D^, and this Dirac operator 
not only specifies the fermion part of the action of the model but also the Yang- 
Mills part Sym(A): there is a definition of Sym(A) in terms of (see ||CC|| and 



references therein). Our discussion above suggests a simple physical interpretation of 
this spectral action principle ||CC| : the logarithmic divergence of the fermion effective 



S YM (A) := const. — S ios (A) (5) 



action is potentially 'dangerous' since it can make the model ambiguous: there is 
no preferred choice for the cut-off, and changing it generates a term proportional 
to Siog (A). However, the fact that S\ og (A) is proportional to the Yang-Mills action 
resolves this problem for the standard (purely vector) Yang-Mills theory on R 4 , as 
discussed above. It therefore is natural to require that the Yang-Mills action is pro- 
portional to the logarithmic divergent part of the fermion effective action in any gauge 
theory models. In particular this suggests the following definition of the generalized 
(vector and chiral) Yang-Mills action in terms of the generalized Dirac operators Da, 

1 

2? 

(for one fermion flavor const. = 24ir 2 ). Eq. (§) shows that for flat Euclidean space 
1R 4 , this definition is equivalent to the one given in ||CC|| . We conjecture that this is 
true for other four dimensional spin manifolds as well. 

The plan of this paper is as follows. We summarize our notation and results in 
Section 2. Section 3 contains a summary of the mathematical prerequisites, i.e., the 
three ingredients of our method mentioned above. The computations of S\ og (A) is 
presented in Section 4 with some computation details deferred to Appendix B. We 
conclude with some remarks in Section 5. Appendix A contains some discussion on 
regularized traces and the noncommutative residue. 

Notation: We write g\ N for the complex N x N matrices and GL^r for the invertible 
matrices in gLy. We sometimes write Iy or / for the identity operator on a vector 
space V but often abuse notation and do not distinguish between cl and c for complex 
numbers. For V, W vector spaces and a an operator on V, we often use the same 
symbols a to also denote the corresponding operator a® Iw an d Iw ® a on V ® W 
and W ® V, respectively. The real part of a complex number c is denoted as 3?c. 



2 Definitions and Results 

For simplicity we assume spacetime M 4 = M 4 with Euclidean signature (the extension 
of our calculation to other four-dimensional spin manifolds should be possible using 
symbol calculus of pseudo-differential operators 0). 
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We consider the Hilbert 

H = L\R i )®Ct pm ®CZ lor (6) 

which has the physical interpretation as space of the 1-particle states of the fermions. 
We also introduce the space T> of functions in 7i which are smooth (i.e., C°°) and L 1 ; 
T> is a convenient dense domain in 7i. 

The Dirac operators of interest to us are of the form 

Da = Y + V v (x) + i-y 5 C v (x)) (7) 

where A = (V,C) (repeated indices u,fi... = 1,2,3,4 are summed over; x = 
(x 1 , x 2 , x 3 , x 4 ) G M 4 ), with d v = and 7^ the Dirac spin matrices acting on Cj pin 
and obeying 

7 V + 7 V = 2rT (8) 
for /i, v — 1, 2, 3, 4, where rf v = i]^ u = diag(l, 1, 1, 1) is the metric tensor, and 

7s : = 7 V7 3 7 4 (9) 

as usual (for the convenience of the reader, explicit formulas for these matrices are 
given in Appendix A.l). 

For simplicity we assume that the functions V u and C v R 4 — > gl^ are regular, 
i.e., they are C°° and vanish like 0(\x\~ A ~ e ), for some e > 0, as \x\ — > 00 (the latter 
condition is to ensure that integrals of regular functions over R 4 absolutely converge). 

In particular, the free Dirac operator is defined by the differential operator 

D = -\rfd v . (10) 

We define the gauge group Q as follows. Let GL^v be the group of all invertible 
matrices in gl^. Then Q is the group of all GL^v- valued functions U on M 4 such that 
U(x) — 1 is a regular function. Note that one can write 

Da = |(1 - 75)7" (-ify + V, + iC M ) + |(1 + 75)7" (-i^ + V, - iC M ) 

where ± iC^ are the chiral components of the gauge field. This representation 
shows that it is natural to consider two kinds of gauge transformations, 

V, ± iC„ - (f/i)- 1 (V^ ± iC M ) c/ ± - Kc/i)- 1 ^^, E4 e . (11) 

For ?7 + = U- = U we denote these as vector gauge transformation, otherwise as chiral 
gauge transformation. 
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Note that Da m Eq. (0) is well-defined on the domain T> C 7i, and we find it 
useful to distinguish this formally self-adjoint differential operator in notation from 
the corresponding self-adjoint extension on 7i which we denote as D^, i.e., (D^/X^) = 
L>Af(x) for all f(x) G V. We also write 

Da = D + A (12) 

where Do is the free Dirac operator (i.e., self-adjoint extension of -Do) and A the 
operator defined by multiplication with the generalized Yang-Mills field 

4 

A(x) = Y,'f(V v (x)+vy fi C v (x)) • (13) 

v=l 

We will compute the fermion effective action S\(A) defined in Eq. ([I]), and we will 
show that it can be expanded as in Eq. (^|). As discussed, Tt\ is a Hilbert space trace 
with an ultraviolet (UV) cutoff A > 0, and Ao is an arbitrary, in general complex, 
parameter makes the argument of the logarithm dimensionless. Moreover, the real 
(positive or negative) parameter m corresponds to a fermion mass and serves as an 
infrared (IR) regulator in our computation. Our main result is an explicit formula 
for Si og (A). 

Proposition: The logarithmic divergent piece S\ og (A) of the logarithm of the (regu- 
larized) determinant of the Dirac operator equals 

Sio g (A) = -L J d*xtr N (±F+ (.F+r + \T, V {T-T - 6m 2 C»C,) (14) 

where tr^ is the usual matrix trace in g\ N and^\ 

3% ■= d » A » - d » A i + '^ A % A ^ A i ■= v » ± c » ( 15 ) 

is the curvature associated with the chiral component A ± of the Yang-Mills field. 
(Proof in Section^ with some details deferred to Appendix B.) 
For C = (no chiral field) we obtain 

S ^ S ( A ) = 7^ [ d^xtx N T^ 

247T Z J M 4 

2 [a, b] := ab — ba 
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with 



•F [iv = d,V u - d v V, + i[V M , V u ], (16) 
which is the standard Yang-Mills action. Note that T^ v = i[D^,D u ] with 

D v := -id, + V v (x) (17) 
the covariant derivative, and similarly, 

^± =i[^±iCV,A,±iCV]. (18) 

It is important to note that for m — 0, S'iog(A) in Eq. ( |14|) this is manifestly 
invariant under all gauge transformations Eq. ([TT]) . For there is also a mass 

term oc C^C^ for the chiral gauge field which is only invariant under vector gauge 
transformations, i.e., only the transformations Eq. ( |TTD with U + = U_ = U . The 
parameter in front of this term is fixed by the fermion mass. There is no similar term 
for the vector gauge field (note that such a term would spoil vector gauge invariance). 

It is interesting to note that the result of our computation in Section [| suggests 
that for manifolds M 4 with boundary <9M 4 , S\ og (A) has an additional contribution 

AS log (A) = -L / d 4 xdnr N J„ (19) 

with 

J, := 2C^[D U , G v \ - 2C V \\D V , C M ] + 2i[D M , CG V \ . (20) 

This is a boundary term (by Stokes's theorem). Note that this term is also invariant 
under vector gauge transformations, and it vanishes if the axial Yang-Mills field C M 
is zero. 

It is also worth noting that, as a by-product, we also obtain the explicit expression 
for the quadratic divergent part of the effective action, 

5(2) (^) = t}—2 I d A xtr N (-V^ + C"C M ) . (21) 
16n 2 J M i 

In contrast to S\ og (A) this term is not gauge invariant (as already mentioned, the 
term oc V^V^ spoils vector gauge invariance)! This highlights the fact that the 
regularization procedure we use it not manifestly gauge invariant but only quasi- 
gauge invariant. It shows that the vector gauge invariance of our result for S\ og 
somewhat remarkable. It is also interesting to note that for = ±C M , S^ 2 \A) = 0. 
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3 Calculation tools 



In this Section we collect the mathematical prerequisites for our computation. We will 
explain the three ingredients for our method: Firstly, a definition of the logarithm 
of operators a in terms of an integral of the resolvent of a. Secondly, a few basic 
definitions for PSDO which imply a simple and elegant formula for the symbol of 
the resolvent of Dirac operators D^. And finally, a definition of a regularized Hilbert 
space trace Tta (corresponding to introducing an UV cutoff A). In the next Section 
we will put these ingredients together and obtain an expansion of the effective action 
as described in the Introduction. 

1. The logarithm of operators. Let a be a bounded operator on a Hilbert space 
H. with norm less then one. Then (I = ly_ is the identity operator) 

log(/ + a) = / - (/ - (/ + sa)- 1 ) , (22) 
Jo s 

as can be seen by a Taylor expansion, 

io g( / + a) = f; = _£;/' £(_)., 

n=l n=l J ° 

interchanging summation and integration, and using the geometric series. 
We take this motivation to define 

where Ao is a some complex number. This representation of the logarithm as integral 
of a resolvent will be convenient for us since there is a simple formula for the resolvent 
of (generalized) Dirac operators, as discussed below. 

2. A. Pseudo— differential operators. Generalities. We summarize some basic 
facts about pseudo-differential operators (PSDO) on M 4 (a discussion for general 
manifolds can be found, e.g., in |H[). We consider PSDO a on TC which can be 
represented by their symbol cr[a](p,x), i.e., a gl 4 <E> gLy-valued functions on phase 
space M 4 x IR 4 defined such that [Q] 

(«/)(*)=/ T^i / d i ye i ^a[a](p,x)f(y) (24) 

J R 4 (27T) 4 J R 4 
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for all f(y) G T> (matrix multiplication is understood; p ■ x = x v p v ). In particular, 
D and A are PSDO with symbols 

a[D ](p,x) =i>:= <fPv, <r[A](p,x) = A(x) . (25) 



Note that Eq. (EJ) implies the following equation which encodes the product of 
operators in terms of their symbols, 

a[ab}(p,x)=[ -p- [ d% e^~y>^a[a](q,x)a[b}(p,y). (26) 

We will encounter PSDO a which allow an asymptotic expansion 

oo 

a[a] ~ ^2a K ^[a\ (27) 

3=0 

where cr K _j [a] (p, x) is homogeneous of degree K—j in p|] and goes to zero like Ip^"- 7 
for \p\ — > oo (\p\ : = y/p ■ p) . We write 

n 

a[a](p,x)=J2^ J [a}(p,x) + 0(\p\ K - n - 1 ) (28) 
for all integers n. Eq. (EB) implies, 

£jj n! 9pi a • • • dp in dx h ■ ■ ■ dx in ' 

This equation allows to determine the asymptotic expansions of a[ab] and cr[a _1 ] from 
the ones of a [a] and a[b}. 

2.B. The symbol of the resolvent. Eq. ([23|) expresses log(D J 4+im) as an integral of 
resolvents of the Dirac operator D^, i.e., of operators [c\I + c^D^) -1 with ci >2 complex 
numbers. We will therefore need the symbol of such a resolvent. To determine this 
we note that 

a[ Cl I + c 2 D A ](p,x) =c 1 + c 2 [p , + A(x)] . (30) 

We then could use Eq. (|29| ) to find the expansion for a[{c\I + C2^>a)~ 1 ]{p, x). We 
now present a useful result summarizing this expansion in a simple formula. 

3 i.e., <jK-j[a](sp,x) = s J (TK-jH(p, x) for all s > and \p\ > 
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Lemma: The following holds for all Ci, c 2 G C, 

(x[(ci/ + c 2 D A ) _1 = (ci + c 2 $ + Da]) a[a](x,p). (31) 

Remark: The proper interpretation of this equation is as follows, 

oo 

<r[(ci/ + c 2 D A ) a](p,x) ~ 2^(-l)" (ci + c 2 #) £>a (ci + c 2 #) <r[a](x,p) 

n=0 

where the differential operators d v in Da = —i^^du + A(x) act to the right on the 
functions A(x) according to the Leibniz rule. We note that we will need this equation 
only for a = I. 

Proof of the Lemma: One can check Eq. ([H|) by using Eqs. (pTD and (|29|), taking 
CiJ+c 2 D^ for a and [ci/ + c 2 D J 4]~ 1 a for b, and inserting Eq. (|30|) . A simpler argument 
avoiding tedious expansions is as follows: Note that by definition, (D^f)(x) = [D + 
A(x)]f(x) for all feV, thus 



(( Cl / + c 2 D A ) af) (x) = ( Cl + c 2 [A) + (o/)(x) = 

d * y jH*-v) ( Cl +C2 y + D() + A( x )]) a [a](p, x)f(y) 



d 4 p 

where we used Eq. (21) and the Leibniz rule. Replacing a in this equation by 



\c,\I + c 2 Da) 1 a, we see that this is equivalent to Eq. fl31"|). (Note that this argu- 
ment implies the interpretation of Eq. fl3l|) as given above!) □ 

Remark: We believe that our expansion in powers of the differential operator Da 
is very natural for at least two reasons. Firstly, since under a vector gauge trans- 
formation, Da — > U~ 1 DaU, such an expansion is close to being manifestly gauge 
invariant (we will discuss this point in more detail below). Secondly, it is natural 
from the point of view of power counting: in contrast to an expansion in A, the n-th 
order term in our expansion includes precisely those local polynomials V n in and 
C v (and derivatives thereof) which all have the same scaling behavior V n — > \ A ~ n V n 
under x — > Ax. 

Remark: Loosely speaking, PSDO are useful since they allow to interpolate between 
Fourier- and position space: generically in quantum theory one deals with operators 
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H on some Hilbert space of L -functions on IR n which are a sum of a free part H 
diagonal in Fourier space|] (H f)(p) = E (p)f(p), and a potential term V diagonal 
in position space, (Vf)(x) = V(x)f(x). The symbol of cr[H](jp,x) is then simply the 
sum of Eo(p) and V(x), which is an attractive feature. The price one has to pay is 
that the symbol of ('nice') functions F of H are somewhat complicated: in a first 
approximation, a[F(H)](p, x) ~ F(E (p)+V(x)) + . . ., but there are correction terms 
. . . depending on derivatives. The Lemma above is a special case of the following 
formula, 

(T[F(H)](p,x)~F(E (p-id) + V(x))l 
nicely summarizing the systematic derivative expansion of functions of H. 

3. Regularized traces and the noncommutative residue. We now define 
the regularized trace which we will use. We first note that due to our technical 
assumptions on the gauge fields all operators a which we will encounter are PSDO 
which have symbols cr[a](p, x) which go at least like (9(|x| _4_e ), some e > 0, for fixed 
p and | a; | — ► oo, and are finite for finite p. Thus 

f d 4 p f 

Tr A (a) : = / 777^7/ d 4 x tr a{a](x,p) (32) 

J\p\<A (27T) 4 J R 4 

where tr is the full matrix trace,(] is well-defined for A < 00, and this defines a 
regularized Hilbert space trace: If a is a trace-class operator then Tr A (a) has a well- 
defined limit A — > 00 which is equal to the Hilbert space trace of a . More generally 
one can consider PSDO a for which Ti\(a) can be expanded as 

Tr A (a) = c W( a )A* + c^-D^A^ 1 + 

+ . . . + c«(a)A + c log (a) log ( ^] + c (0) (a) + 0(A.- 1 ) (33) 



with K some non-negative integer. 



We recall that the noncommutative residue [Wo| of a PSDO a with an asymptotic 



expansion as in Eq. (p7[) can be defined as (see, e.g., Eq. (2.7) in Ref. [| 

1 f d 4 p 



Res(a) :=-! -^-8{\p\-l) I d A x tro"_ 4 [a] (x, p) , (34) 



A f(p) = /k« d n x e zp x f(x) denotes the Fourier transform, 
''including the trace tr„ in gl 4 and the trace trjv in gljv- 
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and for PSDO a as above, 

Res(a) = -ciog(a), (35) 

i.e., the residue is equal, up to a constant, to the logarithmic divergent part of the 
regularized trace of a. (An elementary proof of this latter fact is outlined in Ap- 
pendix A.) 

Remark: In our definition Eq. (|39|) of Tr A we use a sharp cutoff, i.e., 

Tr A (a)= I ^-J(f) I Atm[a](,,p) (36) 

where fit) equals the Heaviside step function 9(1 — t). In principle one could define 
a regularized trace using Eq. ( |56"D and choosing any non-negative, piece-wise smooth, 
function f(t) which vanishes exponentially fast for \t\ — > oo and is such that f(0) = 1. 
For example, the choice f(t) = exp(— t 2 ) would correspond to the standard heat kernel 
regularization. 

We will show in Appendix A that ci og (a) is in fact independent of /. 
Using any such regularization one can define the renormalized trace as the finite 
part of the regularized trace, 

TR(a) := co(ar), (37) 

but this is not quite independent of the regularization: as also discussed in Appendix 
A, changing the regularization function / — > / amounts to changing 

TR(a) -> TR(a) + log(s) q og (a) (38) 

with some constant s > depending on / and /: the logarithmic divergent piece 
accounts for the regularization dependence of the renormalized trace, and this is the 
reason for our interest in it, as discussed in the Introduction. 

Remark: We note Eq. (j3~6D is equivalent to 

Tr A (a) = Tr(P A a), P A := /(|D |/A) (39) 

(using the spectral theorem for self-adjoint operators). This naturally extends the 
definition of Tr A from PSDO to a large class of operators on 7i. More generally, one 
could change the regularization by changing Do — >• D# in the definition of P A , for 
some fixed Yang-Mills field B. One can show that this would change TR(a) by a 
term proportional to Res([log(D B ) — log(D )]a) (see, e.g., Eq. (1.6) in ||CDMP|| ). It 
would be interesting to explore this possibility in more detail. 
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4 Computation of effective fermion action 

In this Section we present the explicit computation of the effective fermion action and 
thus prove the proposition in Section ||. Our computation amounts to a quasi-gauge 
invariant gradient expansion, which is essentially an expansion in powers of the UV 
cutoff A. This allows us to extract, in a simple manner, the quadratic and logarithmic 
divergent pieces which is what we are interested in. 

1. Quasi-gauge covariant expansion. We write 

S A (A)= [ £L [ d 4 *trS(*,p)l (40) 

J\p\<A (A) 

where S(x,p) is obtained by computing the symbol of the operator log(D^ + im)/A — 
log(D + im)/A as explained in the last Section, i.e., 



Six 



p)= [ ^1 ([1 - s + sp + si}- 1 -[1-8 + sf>}- 1 ) = 
Jo s 

ro ° du 

{{l + uf + utfi}- 1 -[l + up]- 1 )- (41) 



we used Eqs. (|22|) and (|31~D, introduced the convenient short-hand notion, 

. p + im —id + A , . 

P ■ = > i ■ = — x , I 42 ) 

Ao Aq 

and changed integration variables, s = u/(l + u). The 1 on the r.h.s. of Eq. ( pE0"D is 
the symbol of the identity operator. As explained in more detail below, S here is to 
be regarded as a differential operators acting on 1. It is straightforward to expand 
the integrand in this equation in powers of ^, 

L 

S = J2(-±) n ~ 1 Sn + n L+1 (43) 

71=1 

where 



and 



s n = (-i)"- 1 [°° ( ui—^— ) (in 

du 1 ( , 1 1 



u 1 + up \ 1 + up; 
13 



is a remainder term. 

In the following we find it convenient to use the short-hand notation and write 

Da = < D l ( 46 ) 

s=0,5 

where 

D° V :=D V D B V :=C V (47) 

and 

Yo--=Y, 7?:=i7"7fe- (48) 

We then define 

where Sj = 0, 5 and i/j = 1, 2, 3, 4. This allows us to write 



J\p\<a {^r Jr* ~~ ~ Jr 4 



(50) 



here and in the following, s is short for (si, . . . , s n ). 

The following Lemma simplifies the computation significantly: it implies that the 
S n for odd integers n all vanish, and that an series expansion in the mass m only has 
non-zero even powers. 

Lemma: The coefficients Ai~. s in Eq. faffi ) are non-zero only for even integers n, 
and they are invariant under m — > —m, i.e., they are independent of the sign of the 
mass. 

(Proof in Appendix B.) 

Remark: We now can explain why we denote our expansion quasi-gauge invariant. 
This is because the operators D s v transform gauge covariantly under a vector gauge 
transformation U, D s v — ► U~ X D S V U . This implies that the differential operators de- 
fined in Eq. (|50D all are gauge invariant. However, the action is a polynomial which is 
obtained by applying these differentiation operators to 1 (cf. Eqs. (40)) using Leibniz 
rule and d v l = 0, e.g., 

D v l = V„(x) 
D Ul C U2 l = -\{d Vl C V2 (x)]+V Ul (x)C U2 (x) 
D U1 D U2 1 = -i[d Ul V U2 (x)] + V Ul (x)V U2 (x) (51) 
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etc. The result will be gauge invariant only if the differential operator in Eq. fl5Tj|) is 
already a polynomial. This happens, e.g., if the differential operators D u only appear 
in combinations [D V ,D^\ and [D u ,C fl ]. This is not obvious. However, we will see 
below that this happens for the terms leading to logarithmic divergent part. 

2. Expansion in powers of the UV cutoff. We now show that our expansion 
above is essentially an expansion in powers of the UV cutoff A. Our computation 
can be simplified by the following argument (this argument is refined and justified in 
detail in Appendix B). As mentioned, m serves as a particular IR cutoff for momen- 
tum integrals. We expect that our result is independent of the precise form the IR 
regularization. Thus we use instead the following, simpler one: we set m = in f> but 
restrict integrations over p to m < |p| < A. We stress that we use this simplification 
in the main text only to ease our presentation, and that it is appropriate only for 
computing the diverging contributions to the regularized determinant: the computa- 
tion of the finite part should be done with the method explained in Appendix B. Below 
we shall see that this simplified procedure gives a IR regularization provided we also 
set Ao = | Ao |/(1 + i0 + ) (a justification of this can be also found in the Appendix B). 
Using then 



1 1 



■[l-^(Ao)- 1 ] 



l + ufi l-u 2 \p\ 2 (A ) 
and rescaling m|p|(A ) _1 — > u{l + i0 + ) we see that M. n in Eq. (f49|) becomes^ 



M^ s = ^- 2 [ A d\p\\p\ 3 ~ n JUp) (52) 

^ J\m\ 



where 



r>00 

J%"» := / duu n ~ l 



n+l 



o vi-Ki + io + )] 2 , 

x (tr„(l - uiKlil - u() ■ ■ - 7 £(1 - u£)) ; (53) 
we used (2n)-*J m ^ A d*pg(p) = (Stt 2 )- 1 j£, d\p\ \p\ 3 (g(\p\0) with 

^■■=^Lw^ m - 1)9{0 (54) 



6 We use the different symbol M n to indicate that these numbers are obtained with a simplified 
IR regularization. 
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the angular average (i.e., integration over the unit sphere in R 4 ). We now see that 
A = |A |/(1 + i0 + ) is needed to specify how to treat the singularity in the it-integral. 
These u-integrals are then finite (see Eqs. (|56| ) and ( |57j ) below). The result we get is 
independent of |A |, as expected. It shows explicitly that our expansion leads to an 
expansion of the action in powers of A. We are interested in A — > oo. In this limit, 
A4 n oc |A| n ~ 4 for n < 4 and cx |m| 4 ~ n for n > 4: the former terms are divergent in 
the UV (i.e., for A — > oo), the latter in the IR (i.e., for m — > 0). It is precisely the 
'boundary case' n = 4 which gives rise to the logarithmic divergence. 

This result is obtained with the simplified IR treatment is correct only in lead- 
ing order in A. In Appendix B we show how to do the computation without this 
simplification, and that 

M n = Mn + O (m 2 A 2 - n ) , n > 2 (55) 

showing that the simplified IR treatment gives the correct result for the diverging 
terms for all n but n = 2. For n = 2 there are corrections oc m 2 log(A/m) which 
contribute to S\ og and which we therefore have to compute exactly. 

3. Computation of diverging parts of the effective action. We now proceed 
to compute the coefficients J n Eq. ( |53D for those terms we are interested in, i.e., for 
n = 1,2,3,4. Using Eq. ( ^2|) this is straightforward: one only needs to evaluate the 
integrals 

/•oo 

Af n , k = / duu^- 1 (1 - Kl + K^)] 2 )^ 1 , (56) 
Jo 

the angular averages (£ Ul • • • £„ fc ), and traces of products of Dirac matrices. The 
integrals in Eq. (||) are (cf., e.g., Eq. 3.251(11.) in gJRJ) 

K,k = -(-l) [n ~ k)/2 B (l±±, ^ + l) (57) 
where B(x,y) = T(x)T(y)/T(x + y). We only need 

^2,0 = ~, A/" 2 , 2 = -- 

M,o = i ^ a = — , ^> 4 4- 

The computation of the traces of Dirac matrices is simplified using the following 
relations 

(1 - u£)Y = 7"(1 + u() - 2uC, ^ = 1,2,3,4 

75^75 = f = K| 2 (58) 
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which follow from Eq. ([]). We also need 

(-0 — 1; (^1^2) — 

) (59) 

and that the angular average for a product of an odd number of components £ Mj is 
zero. Moreover, 

tr7v(7 w 7 M2 ) = V 1M2 
trjv(7 /xl 7 M2 7 /X3 7 At4 ) = 4 (fj^ 1 ^ 2 ^ 3 ^ 4 — ^mi/^^/m _|_ ,^i/m^m2M3^ 

tr^^^VW 4 ) = 4e Ml/ ^ 3/M (60) 

where e^ 1 ^ 2 ^ 3 ^ 4 is the completely antisymmetric symbol with e 1234 = 1. Note that 
tr^v(7^) = tiTv (7^7^7^) = always. 

It is now easy to see that J^ s = J^l^g 3 = 0, thus 

3 1 = S 3 = . (61) 

The simplest non-zero terms are for n — 2. Combining the formulas given above it 
is easy to see that 



1 

^55 = -An = 1, ^50 = ^05 = 0. (62) 



_ 4 „MiM2 

*- y 2;sis 2 ~~ 87T 2 



Thus 



5 2 = A 2 -A / d*x\x N (-D»D li + C^C M ). (63) 

l07T z J 

This is a gauge invariant differential operator. When acting on 1 (cf. Eq. (0)) we 
obtain the quadratic divergent part of the effective action Eq. (ET]) which is not gauge 
invariant. 

As mentioned, is only the leading order contribution to S^. A more careful 
computation without the simplified IR regularization gives (see Appendix B), 

S 2 =S 2 - m 2 log (A") _L f d 4 xtv N (C»C,) + ... (64) 
\ \m\J 8tt z J 

where '. . .' are terms which remain finite for A — > 00. We see that the sub leading term 
which was missed by the naive IR regularization contributes to S\ og . As discussed in 
Section [|, this term is gauge invariant. 
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Table 1: Parameters in Eg. (p5|) where s = (s\, S2, S3, S4). 
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1 
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1 
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S2 
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1 
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S3 
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1 
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S 4 
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1 
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-2 


-2 
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-2 


-2 


-2 





























B s 


-2 
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2 
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-2 


























C s 


2 


-2 





2 


-2 





-2 


2 


























D s 


























— i 


— i 


i 


i 


i 


i 


— i 


— i 



We now turn to the case n = 4 which leads to the logarithmic divergence. All 
relations needed to compute the J^ v ^ v ^ were listed above. The result can be written 
as follows 

jvmwa = I (Ajf 1V2 rf 3Vi + Bsjf^rf^ + C s jf^ri v ^ + Dse UlU2l/3UA ) (65) 
3 

where s = (si, s 2 , s 3 , s 4 ). The numbers Ag, Bg, Cg, As are all given in Table [I]. 
(We have checked this result extensively using the symbolic programming language 
MAPLE.) We note that the numbers Ag, B^, Cg (Dg) all are real (purely imaginary) 
and non-zero only if an even (odd) number of the Sj equal 5. 
Combining these results we find 

^ =i ° s {w\)^L d ' xt "' lVR+v ' 1 (66) 

where 

s 

V I = Y J D^ m D s v \DilD s ^Dil (67) 

s 

with the coefficients given in Table [I|. Vr is a sum of 19 non-zero terms. We now 
claim that it is possible to writeP^ = Vj^x + "Pr,2 where 

V Rjl = -[D^D V ][D„D U ] - [C"\C"][C;,a] + [D^D V ][C„,C V ] 

+[C,C v ][D ll ,D v ] +2[D»,C v ][D li ,C v ] +2[D»,C v ][C ll ,D v ] (68) 
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and 




(69) 




[[D V1 , D U2 ] + [C Ul , C U2 ], [D V3 , C U4 ]] . 



(70) 



(The proof of Eqs. (p5[)-([7D|) are straightforward calculation which we skip.) 



We see that, V RA equals |(^+ {F + Y v + J 7 '^?-)^) with T% defined in Eq. ([15]). 



The remaining terms are linear combinations of commutators! Using the cyclicity of 
the matrix trace we thus obtain 



Remark: Note that Vr and Vi are not differential operators but polynomials (i.e., 
there are no terms (• ■ ^D^). This implies that both these terms are gauge covariant 
which, as we believe, is remarkable. 

5 Conclusions 

The regularization which we used was simple but not manifestly gauge invariant. 
For the result computed in this paper the latter property is irrelevant: since the 
logarithmic divergence is regularization dependent one can compute it using any 
regularization. However, we believe that our method is useful even for computing 
the finite part of the effective action, i.e., S^°'(A) in Eq. (|j). We stress again that 
the simplified IR regularization used in the main text is not appropriate in this 
computation but the formulas given in Appendix B should be used. We conjecture 
that S^(A) computed in this way is gauge invariant. 

As mentioned in the Remark at the end of Section |3], we defined a renormalized 
trace TR|d | using the free Dirac operator D . More general we could use the Dirac 
operator Db with some fixed non-trivial Yang- Mills field B. In particular, we expect 
that the standard (^-function regularization of the logarithm of the determinant of 

should be identical with 



tr N V R , 2 = d"tT N J m tr N Vi = 0. 



(71) 



This implies Eqs. (|TJ])-(^t]) and completes our computation. 



□ 
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were the regularization function is f(t) = exp(— t 2 ). The latter definition has the 
advantage that it is manifestly gauge invariant, but it seems less easy to use for 
explicit computations as ours. It is natural to expect that the difference between the 
latter definition and S^(A) in Eq. (f|) is also proportional to SiogfA). 

Effective action computations are used in many applications of quantum field 
theory. We believe that the methods which we presented should be useful in other 
such contexts as well. 

Acknowledgment: I would like to thank A. Laptev, J. Mickelsson, S. Paycha, F. 
Scheck and K. Wojciechowski for their interest and helpful discussions and S. Paycha 
for comments on the manuscript. This work was supported by the Swedish Natural 
Science Research Council (NFR). 

Appendix A: More on regularized traces 

In this Appendix we outline elementary proofs of some facts about regularized traces 
stated in the main text. 

The logarithmic divergence. We compute the regularized trace in Eq. ( p6|) for an 
operator a with a symbol allowing for an asymptotic expansion as in Eq. fl2"Tp. It is 
easy to see that the contribution of a k [a](p, x) to TrA(a) is 

I™ d\p\\p\ k+ *f(^) I t^j^I-I)/ d*xtra k [a]&x) 
Jo Jri (2vr) 4 J R 4 

where we used the homogeneity of u k [a\. Changing variables, \p\ — > u = \p\/L, and 
comparing with Eq. ([33|) we see that for all k > —3, 

c k+ 4(a) = N k [ ALs(\Z\-l)tT(T k [a](Z,x) (Al) 
jRi (A) 

with N k = J °° duu k+3 f (u) constants depending on /. For k = —4 the computation 
above does not make sense (the constant AL4 diverges), but we can compute c\ og (a) 
as follows. We first subtract from the symbol of a the diverging part which we already 
accounted for and define, 

K+3 

(T± 3 [a](p,x) := a[a\(p,x)-^2a K - j [a](p,x) = (j- A [a\(p,x) + 0(\p\- 5 ). (A2) 

j=0 
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Eq. (^) then suggests that 

dog(o) = lim — f t^j/(¥) / d 4 xtra[a}^ 3 (p,x) . 

A^oo log(A) J R 4 (2tt) 4 J K 4 

Computing this using L'Hospital's rule we obtain 

Piog(o) = lim A f / d 4 xtr<x[a]± 3 (p,x) 

= A lim (7 7^(-/'(t)t) / ^WaUM + e^A- 1 )) ■ 

A ^°° VilR 4 ( 27r ) JR 4 / 

Changing variables etc. as above and using du(—f'(u)) = /(0) = 1 (independent 
of /!) we obtain 

c iog(«)= / T^i^KI- 1 )/ rf 4 xtra_ 4 (a)(^,x). (A3) 

J R 4 (27rJ 4 JlR4 

Recalling Eq. (0) we obtain Eq. fl3J). □ 

Renormalized traces. It is obvious that changing the regularization functions 
f(t) —>■ f(t) = f(t/s) for some fixed s > 0, amounts to changing A — > sA, and thus 
changes — > c^ ^ + log(s) ci og . Thus ([$3|) is obvious for this special case. For more 
general changes f(t) — > /(£) of the regularization function, Eq. (^) can be shown 
using 

/ W^fif) I d**M«]-s(P,z) =Ctog(a) log fAV c (0) (a) + o(A -i) 

Jr 4 (27TJ J R 4 

which follows from our discussion above. 



Appendix B: Computation details 

In this Appendix we present some details concerning our computations discussed in 
the main text. In particular, we give explicit formulas for the Dirac matrices, and we 
also show show how to compute the structure constants M. n in Eq. ( [49|) exactly, i.e., 
without the simplified IR regularization. We also prove the Lemma in Section 4.1 
and Eq. fl55|) , and we give some details about the computation yielding Eq. fl64j) . 
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B.l. Dirac matrices 

A convenient representation for the Dirac matrices is as follows 
a 

where 1 and are the 2x2 unit- and zero matrices and 









M 






o ) 


, 75 = | 


v o 





7 J =( a . Q J )> ^ = 1.2,3, 7 4 =( n ), 75 = ( ^ _^ ) (Bl) 



o n ( o -i \ / i o 



" i = ! i o ; ' a2= \ i o J ' as \ o -i 

the Pauli sigma matrices as usual. 

B.2. Details about the gradient expansion 

We start by rewriting the M. n in a convenient form. We define 

P £ :=I(l + £ i-), e = ± (B2) 



2 V bl 

which are orthogonal projections, P £ P~ £ = and P e 2 = P £ , satisfying P + + P_ = 1. 
We then can write 



which we insert n + 1 times in Eq. (|49f) 



|p|+im 

S = ± * T- "^T~ 



X 



( II 7~ £,-|p|+im ) tll/ {^e 1 l s [Pe2 ' ' ' 7s^£n+l) • 



£i,...,£„+l=± \j=l 1 + « Aq 

We thus obtain 

= E x -- ■ tr * <*Vtf^ • • -C*Vi) (B3) 

ei,...,£„+i=± 

with 

n+l 

In; Elr ., E „ + i = 2Tn;fc , k SUch that ^ Sj = n + 1 - 2k (B4) 
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and 

Rescaling uA/A — > k and introducing £ = |p| / A yields 



(B5) 



Proof of the Lemma in Section 4.1. We note that 

tr (P V^P ■ ■ ■ Y n P ) = ■ T- 

is invariant under Ej —>■ —Ej (since the latter transformation amounts to the variable 
change £ — > — £ in the integral Eq. defining the angular average). Moreover, the 
cyclicity of trace and 75 = 1 implies that T^ s does not change if we replace all P £j 
and 7^ by 75-^.75 and 757^75, respectively Using 7 5 P £ 7 5 = P_ £ and 757*75 = ~l u s 
we obtain T~ s = (— l) n Tr £jS , and using Tz £>s = Tf s this proves that Tf s — and thus 
M. n in Eq. — is non-zero only for even n. 

From Eq. (|B4]) it is obvious that Ej — > — £j corresponds to to — > n + 1 — k, 
and thus Tf s = Tz es implies that we can replace T n . k by [T n;fc + X n;ri+1 _ fc ]/2 in Eq. 
(B3). We can write the latter as a sum of the terms which are even and odd under 



the change the sign of the mass m — ► — m. A simple change of variables shows that 
it-integrals in the odd term 

i \l n . k (m) + J n;n+ i_ fe (m) -X n . k {-m) - 2 n;n+1 _fc(-m)] 
can be written as follows (n even), 

1 f°° -1 Z' 1 \"+ 1 - fe / 1 

4 i-oo dMM " 1 vT+^e+ifiJ v r+^RTTf 

plus the same integral but with k and n + 1 — k interchanged. The latter integrals 
can be computed using Cauchy's theorem: the poles of the integrand are in u = 
— 1/(£ + im/A) and u = l/(£ — im/A) and thus both always in the same half of 
the complex w-plane (upper or lower, depending on the sign of m). Computing 
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the integral by closing the integration path in the half plane where the integrand is 
analytic (which is possible since the integrand vanishes like (9(|m| -2 ) for \u\ — > oo) 
one sees that the integral is zero. This implies M. n (— m) = Ai n (m). □ 



Proof of Eq. (|55|) : Our discussion above implies that we can replace I n - k in Eq. 
O by m n , k = [J n;fe (m) + l n . k (-m))/2. 

We are interested in the terms which diverge for A — > oo. To isolate them it is 
convenient to determine dA4 n /dA. We thus compute 



Ask • 1 

where we introduced the functions 



— m n , k = A 3 - n — 2 i n , k (T) (B6) 



ra+l-fc 

n-1 



Note that the functions I n ,k{v) are well-defined for all real r\ 7^ 0, have a finite limit 
-fn,fc(0 + ) as i] — > 0, and they have series expansions in r/ 2 .[] 

It is easy to see that with the simplified regularization used in the main text we 
can obtain a formula for M^.'s"^ Sn as m Eqs. (|F3^)-(|F34|) but with I n . k replaced by 



^ = A / d\p\\p\ 3 - n I n . tk {0 + ). (B8) 



We thus get 

d / - \ I 



— (j n;fc - X n;fe ) = ^A 3 ~" (/^(f ) - J n;fc (0+)) = ^(m^ 1 -"), (B9) 
which proves Eq. (|55|). □ 



Remark: We now can explain the reason for our choice A = |A |/(1 + i0 + ) in the 
main text: this yields a regularization specifying the otherwise undefined integrals 
-fn ; fc(0), and from Eq. ( |B"9| ) it is clear that this is the regularization yielding a result 
identical with the one obtained with the proper regularization, up to lower order 
terms. 

Computation of £2. For n = 2we need compute .M2 in Eq- (EH) exactly, using the 
formulas given above. 

, n+l-k 



7 To see this note that I n ,k(v) = ®Jo ds ( ^^ j ( s -i+i^ j 
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Similarly as explained in the main text we compute (cf. Eq. (B3)) 



tT„ (Pe^Pe* ■ ■ ■ l? 2 Pe 3 ) = (1 - Sl^-l)* 1 ) . 

Moreover, the integrals defined in Eq. ( |B7| ) for n = 2 and k — 0, 1 are, 

J 2 „(„) = Sft 1 = I _ 3 2 0( 4 N 

,uw; 2(1 + ir]) 2 2 2' w; 



and with Eqs. (jBq) , (|B3|) and (|B4| ) we can compute dA4 2 /dA. Straightforward 
computations yield 

■M^ 2 = <W^ 2 i^2 ( A ^- 2 + ^ log (^|) 41 + O(A )) 

A55 = -Ax) = l, 4? = -2, 4°o ) = 0, (BIO) 
and with Eq. (|50|) we obtain Eqs. (|63f) -( 
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